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(54) Method for measuring modulation parameters of digital quadrature-modulated signal 



(57) The signal to be measured is converted into 
digital form and read into a memory 13 using a trigger 
as the reference, and a signal of a length 1/5 that of the 
signal to be measured, in the vicinity of the trigger, is 
orthogonal-transformed (14). A frequency error of 
the transformed output l 1 ,Q 1 is estimated (15) and is 
used to correct the output l|.Q| (16). The corrected out- 
put l 2 .Q2 is subjected to processing by a receiving filter 
(17). and the filtered output l 3 ,Q 3 is synchronized with a 
PN code to obtain an idea] signal R ( ,Rq and a phase dif- 
ference x 1 is calculated. The signal l 3> Q 3 and the ideal 
signal R|.R Q are used to obtain a frequency error 02 
and a phase error t 2 of the signal 13,03. All signals to be 
measured, stored in the memory 13, are subjected to 
orthogonal transformation (23) so that the frequency 
error a^+c^ and the phase error t 1 +t 2 are removed, 
and the transformed output l 5 ,Q 5 and the signal R|,Rq 
are used to calculate modulation parameters as is the 
case with the prior art. 
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Description 

BACKGROUND QF THE INVENTION 

5 The present invention relates to a method for measuring modulation parameters of the transmitted signal quadra- 

ture-modulated by a preknown code (data), such as a carrier frequency error, the modulation accuracy, an IQ origin off- 
set an amplitude error, a phase error, the waveform quality and a time alignment error, as is the case with the 
transmitted signal from a base station defined by, for example, the standard of CDMA (Code division Multiple Access) 
digital mobile communications (EIA/TIA/IS-95 Mobile Station -Base Station Compatibility Standard for Dual-Mode 

10 Wideband Spread Spectrum Cellular System, which will hereinafter be referred to as IS-95 standard). 

The output signal from a conventional transmitter of the QPSK modulation system needs to be demodulated accu- 
rately on a bitwise basis. On this account, criterion for assessment by the standard are rigid and the output signal from 
the actual transmitter to be measured is also a relatively clean signal of a waveform not seriously distorted or disturbed, 
that is, an easy-to-demodulate signal. Since the signal of the system of the IS-95 standard is a spread spectrum signal, 

75 however, it is not required to be demodulated accurately bit by bit. That is. the QPSK signal presents no problem even 
if it is a signal of a distorted or disturbed waveform such that a bit error arises when it is demodulated after quadrature 
detection. Hence, it is necessary that a parameter measuring method for determining the waveform quality of the trans- 
mitted signal from the base station, which meets the IS-95 standard specifications, be able to determine the waveform 
quality even if the waveform is distorted or disturbed more than in the case of the QPSK modulated signal. 

20 In view of the above, p is defined by the IS-95 standard as a waveform quality parameter, not as the modulation 
accuracy. The parameter p is obtained by normalizing the value of correlation between the modulated high-frequency 
signal from the transmitter and an ideal modulated signal by their powers. When the both signal match, the parameter 
p assumes the maximum value 1 . 

In the CDMA system the coincidence of timing among base stations in the communication system is indispensable 

25 because of the property of the spread spectrum signal. As regards the transmitted signal from the base station, it is nec- 
essary that the timing for sending the high-frequency signal coincide with certain timing of the standard time in that base 
station. The tolerance of this timing is also standardized by the IS-95 standard. This timing is defined by the difference 
between a trigger that is sent from the base station and specific timing of a PN pattern of the high-frequency signal 
transmitted from that base station. This trigger is sent at the beginning of the PN pattern period in the standard time in 

so the base station. This timing error is called a time alignment error. 

The CDMA signal is spread by a high spreading chip rate rather than by the transmission bit rate. Hence, in the 
assessment of a modulated signal of a duration corresponding to transmission bits of the same number as that of the 
bits subject to assessment in the conventional standard for measuring the modulation accuracy of the QPSK modulated 
signal, the number of chips will be appreciably large even if the duration is the same. In the conventional scheme for 

35 measuring the modulation accuracy of the QPSK modulated signal using the chip number as the symbol number, the 
number of symbols to be measured is large and the computing time is long accordingly. By the way, in RDRSTD-28 that 
is the PHS standard in our country, a signal of a 1 -burst-600-^iS period is subject to assessment and the symbol rate in 
this case is around 200 kHz, whereas in IS-95 the chip rate is approximately 1 .2 MHz. If a signal of the 600-nS period 
is subject to assessment as in the case of the PHS standard, the number of symbols to be estimated will become about 

40 six-fold. 

In the measurement of the waveform quality, it is standardized that the transmitted signal to be measured is applied 
to a complementary filter to eliminate intersymbol interference that is caused by a filter and a transmitting phase equal- 
izer at the transmitting side. 

To this end, the complementary filter is defined to provide a Nyquist filter characteristic when the sending filter, the 
45 transmitting phase equalizer and the complementary filter are connected in cascade (connected in series). 

in the CDMA system the use of such a complementary filter characteristic in a filter at the receiving side provides 
the optimum receiving characteristic. Consequently, this complementary filter serves as a receiving filter of the meas- 
uring apparatus. 

The complementary filter serves as a narrow band-pass filter for the high-frequency signal but a low-pass filter for 
so the base band signal. In the measurement of the waveform quality, processing by this filter is carried out in accordance 
with the carrier frequency of the input signal, but in this instance, if the carrier frequency is not accurate, there is a fear 
of the input signal being partly cut out because the filter is narrow-band. 

Further, the receiving filter is equipped with group delay frequency characteristics that cancel those of the transmit- 
ted signal. On this account, it is necessary in the measurement of the waveform quality, too, to measure the timing 
55 (phase) of the input signal after the filter processing. In this case, a major or serious error arises in the timing measure- 
ment unless the carrier frequency is accurately estimated. 

Conventionally, this filter processing is repeated a number of times to repeat the estimation of the carrier frequency 
as disclosed in the specification of U.S. Patent No. 5,187,719 (issued february 16. 1993), especially in Fig. 15 and the 
associated disclosure. The timing (phase) of the input signal is estimated only once prior to the filter processing since 
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the traditional receiving filter has flat group delay frequency characteristics. Thus, the waveform quality measuring 
method using the receiving filter is also required to reduce the computing time. 

In the measurement of the waveform quality standardized by the IS-95 standard, it is necessary to synchronize the 
spreading code of the measuring apparatus with the spreading code of the input spread spectrum signal, 
s One of channels for the transmitted signal from the base station of the IS-95 standard is called a pilot channel. The 

signal of the pilot channel is a signal QPSK-modulated by PN codes of in-phase and quadrature components (I) and 
(Q). The PN codes of the in-phase and quadrature components (I) and (Q) are called gold codes and their cross corre- 
lation is appreciably low. 

To synchronize the spreading code of the measuring apparatus with the spreading code of the pilot signal, it is cus- 

io tomary in the prior art to convert or transform the input pilot signal into a sequence of digital values, with the sampling 
frequency set eight times higher than the chip frequency of the spreading code (PN code) of the pilot signal, i.e. a quad- 
ruple oversampling. Since the spreading codes are synchronized by calculating the correlation between them while 
shifting the spreading code generator sample by sample, the correlating operation is appreciably time-consuming. 
In the measurement of each of the afore-mentioned parameters, its symbol timing must be estimated. Convention- 

15 ally, the symbol timing is estimated using an M-phase (M being an integer equal to or greater than 2) PSK complex base 
band signal derived from the input spread spectrum signal and demodulated data (phase data) demodulated from the 
complex base band signal. This will be described below with reference to Fig. 1 . 

For example, a QPSK complex base band signal r(t) is input via an input terminal 1 1 1 into a sampling circuit 1 12 
and a differentiation circuit 113. In the sampling circuit 1 12 the signal r(t) is sampled by the oscillation output from a voft- 

so age-controlled clock (VCC) generator 114 and the output from the differentiation circuit 1 14 is also sampled by the out- 
put from the voltage-controlled clock generator 1 14 in a sampling circuit 120. The sampled outputs from the sampling 
circuits 112 and 120 are fed to multipliers 115 and 116, wherein they are multiplied by an ideal signal (a reference sig- 
nal) exp(-jek) (Ok being a k-th phase of demodulated data) from a signal generator 117. These multiplied outputs, that 
is, deviations from ideal data, are added together by cumulative adders (integrators) 1 18 and 119; namely, they are 

25 averaged. The outputs from the cumulative adders 1 18 and 1 19 are multiplied by a multiplier 121 to detect a deviation 
from ideal symbol timing and its real part is detected by a circuit 122. The detected output is used to control the clock 
signal from the voltage-control led clock (VCC) generator 1 1 4, putting the output from the VCC 1 1 4 into synchronization 
with the symbol timing. 

A brief description will be given of the estimation of the symbol timing of the QPSK complex base band signal by 
30 the maximum likelihood estimation method heretofore employed although this is disclosed in "Digital Communications," 
ProaWs McGraw-Hill. 

In the first place, a logarithmic likelihood function Al (<|>.t) can be set as follows: 



35 



A l (()>,t) - Re[exp(-jWfr(t)R*(t-c)dt] 

where r(t) is the received signal (complex base band signal). R(t) is a reference signal, -c is a time lag, <fr is the initial 
phase of the carrier, T is the measuring time, f is an integration from 0 to T and * is a complex conjugate. In the maxi- 
mum likelihood method, the time lag x is calculated so that the logarithmic likelihood function {$,x) becomes maxi- 
mum. 

40 That is. b that satisfies the following equations is eliminated and the equations are solved on t. 

dAjd$=0 dA Jdx = 0 
From the above equations t is calculated to satisfy the following equation. 

Re[Z(t) • (dZ*(x)/dx)] = 0 0) 

where 

Z(-0 = Mt)R*(t-x)dt (2) 
where J indicates an integration from 0 to T and 

R*(t-0 = ig(t-T-kT c )exp(-je k ) ( 3 ) 



45 



50 



55 



In the above, g(t) is an impulse response characteristic of a Nyquist filter (g(t)=0 for |MT s |>0), e k is the k-th phase of 
demodulated data, T c is the chip interval, T s is the sampling period and £ is the summation from k=0 to a value corre- 
sponding to T. Substituting them into Eq. (1) gives 
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Re[Z(x) • dZ*(T)/ax] = Re[X{exp(-j0k) • Y k (x)} • S{exp(-jek)c?Y k (i)/dz}] = 0 (4) 

where 

5 Y k (t)»fr(t)g(t^-kT c )dt (5) 

In the above, / is an integration from 0 to T. 

That is, Z{exp(-j6k) • Y k (x)} corresponds to the output from the circuit 118 in Fig. 1 and X{exp(-jek)aY k (T)/ac} corre- 
sponds to the output from the circuit 1 1 9. The VCC 1 1 4 is controlled so that the real part of the products of these out- 
10 puts is reduced to zero, and the sampling timing of each of the sampling circuits 112 and 120 is estimated from the 
symbol timing of the input QPSK complex base band signal. 

The configuration of Fig. 1 is to estimate the symbol timing by the maximum likelihood estimation method. 

The Fig. 1 configuration is formed by hardware, and in the case of implementing it by software, no highly accurate 
solution could be obtained because the input signal is a discrete time signal. To obtain the solution with high accuracy, 
15 the sampling rate (the sampling speed) must be increased-this inevitably increases the amount of processing by an 
interpolation filter and hence increases the processing time. 

In the estimation of the carrier frequency in the conventional parameter measurements, in the case of a four-phase 
(n=4) PSK signal, that is, in the case of the QPSK signal, the intermediate-frequency (or high-frequency) QPSK signal 
in the digital form is raised to the fourth power to thereby remove the modulated signal. That is, the QPSK signal 
20 assumes any one of four phase points sequentially displaced nJ2 apart in accordance with the modulated signal, so that 
by raising the QPSK signal to the fourth power, any phase states become integral multiples of 2k and the modulated 
signal is removed. 

The signal thus raised to the fourth power is subjected to fast Fourier transform, then a frequency component that 
provides the maximum peak is extracted from the result of the Fourier transform and frequency divided down to 1/4, and 
25 the frequency of this divided output is estimated as the carrier frequency of the input QPSK signal. This estimated fre- 
quency is used to convert or transform the input QPSK signal to a complex base band signal through Hilbert transfor- 
mation or orthogonal detection, besides the estimated frequency is used to estimate the deviation of the QPSK signal 
from the standard carrier frequency. 

The input QPSK signal is a digital signal, and if its sampling frequency f 3 does not satisfy a condition 2f m <f s with 
30 respect to the maximum frequency f m of the QPSK signal, aliasing occurs and original information cannot be retained. 

Since the QPSK signal is raised to the fourth power, its maximum frequency is four-fold, that is, 4f m . If the sampling 
theorem is not satisfied in this instance, aliasing occurs and the carrier frequency cannot accurately be estimated. 

It is therefore an object of the present invention to provide a method with which it is possible to correctly measure 
modulation parameters of the input signal even if it is prone to produce a demodulation error. 
35 Another object of the present invention is to provide a method which permits accurate measurement of modulation 
parameters of the input signal with a small computational complexity and in a short time. 

Another object of the present invention is to provide a method which enables modulation parameters of the input 
signal to be measured by a signal of a shorter duration. 

Another object of the present invention is to provide a modulation parameter measuring method which permits 
40 rapid synchronization of a spreading code with the input signal. 

Another object of the present invention is to provide a modulation parameter measuring method which enables 
symbol timing to be estimated by software in a short time. 

Still another object of the present invention is to provide a modulation parameter measuring method which permits 
an accurate estimation of the carrier frequency. 

45 

SUMMARY OF THE IN VENTION 

According to the present invention, a input digital quadrature-modulated signal is subjected to orthogonal transfor- 
mation to a first complex base band signal In a first step. In a second step a frequency error, an initial phase error and 

so a timing error of the first complex base band signal are estimated and an ideal reference signal is calculated which is 
synchronized with the first complex base band signal. In a third step the input digital quadrature-modulated signal is 
subjected to orthogonal transformation to a second complex base band signal so that the frequency error, the initial 
phase error and the timing error are corrected. In a fourth step a frequency error and an initial phase error of the second 
complex base band signal is estimated from the second complex base band signal and the ideal reference signal. In a 

55 fifth step the estimated frequency error and the initial phase error of the second complex base band signal is corrected 
to provide a third complex base band signal. In a sixth step a modulation parameter is estimated from the third complex 
base band signal and the ideal reference signal. 

In a seventh step the frequency error of the first complex base band signal estimated in the second step is cor- 
rected to provide a fourth complex base band signal, which is subjected to processing by a receiving filter or comple- 
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mentary fOter to obtain a fifth complex base band signal. In an eighth step the above-mentioned ideal reference signal 
is generated which has its symbol synchronized with that of the fifth complex base band signal, and an initial phase 
error of the first complex base band signal is estimated from the phase difference between the ideal reference signal 
and the first complex base band signal. 
5 In a ninth step a frequency error and an initial phase error of the fifth complex base band signal are estimated from 

the fifth complex base band signal and the ideal reference signal. The sum of the frequency errors estimated in the sev- 
enth and ninth steps is the frequency error estimated in the second step, and the sum of the initial phase errors esti- 
mated in the eighth and ninth steps is the initial phase error estimated in the second step. 

In the second step, symbol decision points of quadrature modulation are estimated from the first complex base 
10 band signal sequence, then samples near the estimated symbol decision points are extracted from the first complex 
base band signal sequence, and the correlation between the extracted samples and spreading codes is calculated. In 
this instance, the correlation calculation is conducted only for the real part (in-phase component) or imaginary part 
(quadrature component) in the first complex base band signal sequence. 

In the estimation of the symbol decision points, the cross correlation between the first complex base band signal 
15 and demodulated data is calculated, then the cross correlation is subjected to filter processing using three Nyquist filter 
characteristics approximated as functions of symbol timing x, respectively, and the results of these filter processing are 
used to calculate the symbol timing. 

In the second step, the first complex base band signal is raised to the n-th power, then the signal raised to the n-th 
power is subjected to discrete Fourier transform, and the maximum peak component is extracted from the result of 
20 transform. The frequency of the maximum peak component is divided down to 1/n, and the divided frequency and the 
standard carrier frequency of the input n-phase PSK signal are added together to estimate the carrier frequency of the 
input signal, n is an integer equal to or greater thgn 2 and usually satisfies n=2 m (where m=1, 2, ...). 

BRIEF DESCRIPTION OF THE DRAWINGS 

25 

Fig. 1 is a block diagram showing the functional configuration for a conventional symbol timing estimation; 
Fig. 2 is a block diagram illustrating the functional configuration of a measuring apparatus employing the method 
for measuring modulation parameters of a digital quadrature-modulated signal according to the present invention; 
Fig. 3 is a block diagram showing the functional configuration of a frequency estimating section 15 in Fig. 2; 
30 Fig. 4 is a block diagram showing the functional configuration of a PN synchronized ideal signal generating section 
18 in Fig. 2: 

Fig. 5 is a diagram for explaining the operation of the section 18 shown in Fig. 4; and 

Fig. 6 is a block diagram illustrating the functional configuration of a signal delay estimating section 22 in Fig. 2. 

i 

35 DESCRIPTION OF THE PREFERRED EM BODIMENT 

In Fig. 2 there is illustrated the functional configuration of the measuring apparatus according to an embodiment of 
the present invention. This embodiment shows an application of the invention to the measurement of quality parameters 
of a pilot signal that is transmitted from a base station in the CDMA digital cellular telephone system. The I component 

40 of the pilot signal is a QPSK-moduiated version of the carrier by certain PN pattern data and the Q component a QPSK- 
modulated version of the carrier by another PN pattern data. The pair of such PN patterns is called a gold code, and 
the cross correlation between the two PN patterns is uniformly low without scattering. 

The pilot signal (hereinafter referred to as a BS pilot signal) transmitted from the aforementioned CDMA base sta- 
tion is converted to an intermediate-frequency signal, which is input via an input terminal 11 into an AD converter 12, 

45 wherein it is converted to a digital signal, and the digital signal of the period to be measured is once stored in a memory 
1 3. The sampling frequency in the AD converter 1 2 is set to be eight times higher than the chip frequency of a spreading 
PN pattern of the BS pilot signal. The memory 13 is a ring buffer, which is designed so that when a write or read reaches 
the last address, it starts again at the address 0, and the write in the memory, an address increment and write halting 
are controlled by a trigger signal that is fed from the outside. That is. the digital signals that are written into the memory 

so 13 are before and after the trigger signal. Some of the stored contents of the memory 1 3. for example, signals of at least 
2048 samples (in terms of the chip number, 1/5, i.e. about 400 samples) are read out before and after the trigger signal. 
In this instance, the number of samples of the digital signal is significantly smaller than the data length of the period for 
measurement defined by the standard and it is not clear whether the data is that in the period, but it is true that the data 
is that neighboring the beginning of the PN period. The signal thus read out of the memory 13 is fed to an orthogonal 

55 transform section 14. wherein it is subjected to orthogonal transform (orthogonal detection) to a first complex base band 
signal l 1a Q n through multiplication by sin cat and coscot. 

The first complex base band signal I, is fed from the orthogonal transform section 1 4 to a frequency estimating 
section 15, wherein a carrier frequency error of the input signal is estimated. The input signal is a QPSK signal that 
assumes four phase states, and its carrier frequency is obtained by such a method as shown in Fig. 3. That is, the first 
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complex base band signal is raised to the fourth power in a fourth power calculating section 124 to degenerate the mod- 
ulated component to obtain one phase state, then the signal raised to the fourth power is subjected to fast Fourier trans- 
form by an FFT section 125, and that frequency in the FFT output at which the maximum peak is obtained is frequency 
divided down to 1/4 in a frequency dividing section 126 and the standard frequency of the CDMA digital cellular teie- 

5 phone system is added to the frequency-divided output in an adder section 127 to obtain the carrier frequency. The 
sampling rate of the first complex base band signal is 4.9152 MHz, four times higher than the chip rate, the number of 
FFT points in the FFT section 125 is 1024, frequency resolution is 4.8 kHz and the frequency estimation accuracy is 
±2.4 kHz. The output frequency of the frequency dividing section 126 represents a deviation from the carrier frequency 
defined by the standard, i.e. an estimated frequency error. Since the data used for the estimation of the frequency error 

to is data neighboring the beginning of the PN period, even if it is displaced from the period to be actually measured, the 
carrier frequency or timing error does not change in a short time in the transmitter, and hence it is regarded as the car- 
rier frequency error in the measuring period. 

The thus estimated frequency error is used to correct the first complex base band signal I^Qj in a frequency 
correcting section 16 (Fig. 2). This correction is made by the following calculation. 

15 

l 2 +jQ 2 =(i 1 +jQ 1 )exp(-ja> 1 t) 

The thus corrected complex base band signal l 2 ,Q 2 is subjected to delay complementary filtering in a receiving filter 
17. That is, in the base station the signal to be transmitted therefrom is subjected to filter processng so that it has a 
20 characteristic which will cancel the group delay frequency characteristics in the mobile station to which the signal is to 
be transmitted. For this reason, the receiving filter 17 cancels the group delay frequency characteristics of the filter in 
, the base station. 

The filtered complex base band signal I3.Q3 is provided to a PN synchronized ideal signal generator section 18, 
wherein it is synchronized with a locally generated PN pattern. Prior to the synchronization, symbol decision points are 

25 estimated. The synchronization of the entire sampling data (in this example, the number of samples is eight times as 
large as that of one chip) with the PN pattern data inevitably increases the computing time. To avoid this, in-phase and 
quadrature PN codes are used to pre-estimate symbol decision points of the QPSK modulation and only sample points 
closest to the symbol decision points, that is. data every eight samples, are synchronized with the PN pattern data. The 
symbol decision points are estimated as described below. As shown in Fig. 4, the in-phase and quadrature components 

30 of the corrected complex base band signal l3(k),Q 3 (k) is squared in multiplier sections 132 and 133, respectively, and 
the squared outputs are added together in an adder section 134 to obtain an amplitude squared sequence of the com- 
plex base band signal l 3 (k), Q3(k), which is subjected to discrete Fourier transform in a discrete Fourier transform sec- 
tion 135. A symbol (chip) frequency component in the transformed output is extracted and the arctangent of the 
extracted component is calculated in an arctangent section 136. From the arctangent output a sequence of zeroes is 

35 extracted, and the zeroes are alternately separated into two sequences, the one of which is estimated as a sequence 
of symbol decision points. 

In the above example the sampling frequency of the complex base band signal sequence is chosen to be eight 
times higher than the chip frequency. In this instance, a sequence of l 3 (k) 2 +Q 3 (k) 2 is multiplied by eight cosine values, 
for example. 1. 0.707, 0. -0.707. -1. -0.707. 0 and 0.707 in a sequential order and a sequence of ^(k^+Q^) 2 is simi- 

40 larly multiplied by eight sine values one after another. By this, the chip frequency component is obtained which is the 
result of the discrete Fourier transform. 

The samples near the chip decision points estimated in a symbol decision estimating section 131 are decimated in 
a decimating section 137 from the real part (the in-phase component) l 3 (k) in the complex base band signal sequence. 
Assuming that samples in the sequence of the real part l 3 (k) are at the positions indicated by crosses X in Fig. 5A and 

45 the estimated symbol decision points at the positions indicated by white circles O* that one of the samples at the posi- 
tions X which are closest to each position O is taken oat as shown in Fig. 5B. In the case where the sampling fre- 
quency is chosen to be an integral multiple of the symbol frequency, once the difference between a first symbol decision 
point and the sample point corresponding thereto is determined, it is enough to take out the subsequent samples at 
intervals of 4T S (where T 3 is the sampling period). The correlations between a sample sequence l 3 (4k) of the real part 

so l 3 (k) near the first symbol decision point and spreading codes R f (4k) and R Q (4k) from a spreading code generating sec- 
tion 140 are calculated in correlation calculating sections 138 and 139. The results of the calculations are squared in 
multiplier sections 141 and 142, and the squared outputs are added together in an adder section 143. As the result of 
this, a correlation value A'(4k) between the real part ^(k) and the spreading codes R,(4k) and Ro(4k) is provided from 
the adder section 143. The phases of the spreading codes that are generated in the spreading code generating section 

55 140 are shifted chip by chip so that the correlation value A'(4k) becomes maximum. When the correlation value A'(4k) 
becomes maximum, the spreadng codes R ( (4k) and R Q (4k) are synchronized with the input l 3 (4k). Incidentally, the cal- 
culation of the correlation value A'(4k) needs only to be conducted over ±100 chips of the real part l 3 (k). 

In the above-mentioned pilot signal, the correlation between the in-phase and quadrature PN codes is uniformly 
low, so that even if the phase of the complex base band signal is not correct, correct synchronization could be accom- 
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pfished by synchronizing only the in-phase PN code with the complex base band signal. Hence, the spreading codes 
R,(4k) and Ro(4k) may also be similarly synchronized with the imaginary part Q 3 (4k) in place of the real part l 3 (4k). 

In a different case, for example, in the case of synchronizing the PN codes with a spread spectrum signal obtained 
by spreading the QPSK modulated signal by a spreading code, it is necessary only to calculate complex correlations 

5 between l 3 (4k), Q 3 (4k) and R|(4k), Ro(4k). In this instance, too, since the calculation for correlation is conducted at 
intervals of the symbol sequence (4T S in the above example), the computational quantity for calculating the correlation 
is smaller than in the case of calculating at intervals of the sample period. This estimation can be made with accuracy 
in the symbol decision estimating section 131 by inputting again the data l 3 (4k) and Q 3 (4k) near the symbol decision 
points once estimated into the symbol decision estimating section 131 to estimate the symbol decision points in the 

10 same manner as described above. 

By the PN synchronization described above, the ideal signals R|,Rq can be generated. 

Turning back to Fig. 2. the ideal signal R|.Rq thus produced in the PN synchronized ideal signal generating section 
1 8 and the measured signal (complex base band signal subjected to filter processing) I3.Q3 is compared in a frequency 
and phase error estimating section 19 to estimate a carrier frequency error and an initial phase of the input signal and 

15 the phase difference between the ideal and measured signals at each of their corresponding symbol decision points is 
also calculated. If the phase difference exists, it linearly increases with the lapse of time. That is, representing time on 
the abscissa and the phase difference on the ordinate, the phase difference is a linear function of time and its intercept 
or segment, that is, the phase difference at time 0 is the initial phase and the inclination to the time axis is the carrier 
frequency error. Letting the sample number and the phase difference be represented by x and y, respectively, the incli- 

20 nation b of a phase difference line by the least squares method is given by the following equation. 

b = (Zyj-lXj -nSx, ♦y i )/(Xx i - Sx, -nSXj 2 ) 

where S indicates the summation from i=1 to n and Xj is the sample number (time). Hence, representing Xj by i, the 
25 above equation becomes as follows: 

b = (Syj • Si - nSi • y ( )/(Si • Si - nSi 2 ) 

where 2 indicates the summation from i=1 to n. Si m n(n+1)/2 , and Si 2 - n(n+1)(2N+1)/6 . 

30 According to the definition of the carrier frequency calculation, a positive frequency error exists when the phase 
rotation is clockwise on the IQ plane, and since the direction of the sign of the inclination b is positive, the Sign is 
reversed. The inclination b is regarded as a carrier frequency error 002. The initial phase is made more accurate by 
deriving it from the above-mentioned intercept or segment, but the initial value of the phase difference is used because 
the estimation needed in this case may be rough. 

as The initial phase e 0 and the carrier frequency error a>2 thus obtained are used to correct the filtered base band sig- 
nal I 3 ,Q 3 by the following equation in a frequency and phase correcting section 21 . 

I 4 + jQ 4 = (l 3 + jQ 3 ) exp (|© 2 t + e o) 

40 Next, the amount of delay of signed is estimated in a signal delay estimating section 22. The amount of delay of the 
PN pattern behind the trigger is already known and the sampling point closest to the symbol decision point is also 
known. The signal delay estimating section 22 estimates the signal delay t 2 in such a manner as described below. The 
delay %z can be estimated by calculating the parameter clock phase x t the carrier phase 4> and the carrier frequency f 
so that the logarithmic likelihood function given by the following equations is maximized. 

45 

A(<|». f. t) = Const. {exp(-j<t>)C(f. x) + exp(j<»C*(f, t)} (1) 

C(f, t) = f TO 0 Z(t)exp(-2j7tft)R*(t-x)dt (2) 

so where Z(t) is the measured signal. R(t) is the reference signal and T Q is the measuring time for estimating the parame- 
ters. 

The continuous signal R(t) and the discrete signal R(k) bear the following relationship. 

R(t) = R(kT s ) 

55 

where T 3 is the sampling period. 

The principle of estimation similar to that by Eqs. (1) and (2) is disclosed, for example, in Proaks, "Digital Commu- 
nication; 2nd ed.. pp. 333, Eqs. (4, 5. 71), McGraw-Hill. 1989. More specifically, the following simultaneous equations, 
in which equations differentiated by the respective parameters are set at zeros, are solved for the parameters 4>, f and x. 
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exp(-j<»C(f, x) • exp(j<»C*(f, x) = 0 (3) 

exp(-j4»)dC/dT + expQWaCVar « 0 , W 

axp(-j*)dC/af + expfl+)dC*/af--0 ( 5 ) 
By simultaneously solving these equations, they can be modified to the fallowing simultaneous equations without ft. 

a|C(t, x)| 2 /&x =0 (6) 

d|C(f. t)| z /af = o (7) 

From Eq. (6) x is calculated setting f»0. Concrete steps of the above calculations will be described below. 
In the first place, C(0. t) is expressed with the integration replaced by the summation as follows: 

C(0. x) = T 8 22(kT 9 ) - R*(kT 8 -r) (8) 

where £ indicates the summation from k=0 to K-1 and T a =T c ;8, T c being the chip period. The measuring time 
• T 0 =KT g . The continuous signal R(t) is expressed by the phase e n of the ninth chip as follows: 

R(t) » SexpQ9 n )-s(t-nTc) (9) 



where Z indicates the summation from n=-oo to oo and s(t) is the characteristic of the base band filter, in practice, the 
summation of Eq. (9) must be effected in a finite or limited period. To this end, the base band filter is assumed to be 
25 zero-phase and to have a response time T F That is, 

s(t)=0(|t|>T F /2) (10) 
In this case, the ideal signal can be expressed as follows: 

R(kT 8 - x) = I n expG6 n )s({8n-k}T s + a) (11) 

Then, using Eq. (1 1). Eq. (8) can be expressed by the following equation by the use of the discrete signals Z(ty and 
R(k).' 



30 



35 



M 

C(0.t)=T s £ £Z({8n-m}T 3 )exp(-j9n)s(mT 3 +x) (12) 



40 Letting T F =(2M+1 )T a , and the measuring time be represented by KT S (k*=0 to K-1 ) , the 6 n must be measured for a time 
n=- M/8 to (K+M-1J/8 . Next, s m (x)-s(mT s +x) is approximated by a quadratic expression of x as follows: 

s m (x) = a m + b m T + c rn x 2 (13) 

4$ However, since x to be estimated is within the range of |x|<T 3 /2 , the approximation needs only to be accomplished 
within this range. Using this, C(0, x) is given by the following equation. 

M 

C(0.x)=T 9 £ £ Z « 8n " m J T J ex PH e nXam-*m^+Cm^ J (14) 
SO m— M n 

= T s (A+B-c + Cx 2 ) (15) 
ss A, B and C in the above are given by the following equations. 

A=££Z({8n-m}T 8 )exp(-je n )a m (1 6) 
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B^ZZpmlTJexpl-jejb,, (17) 



C=X2Z({8n-m}T s )exp(-je n )C rn ( 1{J ) 
m n 

Substitution of Eq. (15) into Eq. (6) gives the following equation for x. 

Re[C(0, x)aC*(0, x)/dx] = (T 3 ) 2 Re[(A + Bx + Cx^XB' + 2C*x)] =0 (19) 

This is a cubic equation for x, but H a linear approximation is used on the assumption that x is small, an equation for 
calculating the clock phase x is given as follows: 

15 x*-Re[ABT/(|B| 2 + 2Re[AC*]) (20) 

This x is the aforementioned x 2 . 

That is, as shown in Fig. 6, the cross correlation between the ideal reference signal R|(8k)+jRQ(8k) and the meas- 

20 ured signal U(k)+jQ 4 (k) is calculated in a cross correlation calculating section 151. then the calculated result is filtered 
by filters 152, 153 and 154 of filter coefficients a m . fc^ and c m , respectively, to conduct the calculations of Eqs. (16), (1 7) 
and (18) to obtain A, B and C, which are input into a timing error calculating section 155 to calculate Eq. (20) to obtain 
the timing error x 2 . The filter coefficients a^. b m and c m are calculated, for example, as described below. To prevent 
intersymbol interference in the Nyquist filter, the time from the peak of its impulse response to the first zero point needs 

25 to be T c (the symbol period, in this example, the chip perio d).Letting a roll-off coefficient be represented by a. the 
impulse response of this filter is given as follows: 



sinfrVTJ cos(Ttat/T c ) 

S(t)= (Ktrrj * i-(2atn- c 

30 

This impulse response is approximated by Eq. (13) to obtain the filter coefficients a m , b m and c m at that time. 

In an orthogonal transform, frequency, phase, x correction and receiving filter section 23, the parameters estimated 
so far. that is, the carrier frequency error, the initial phase and the signal delay x, are used to perform orthogonal trans- 

35 formation of all the intermediate-frequency digital signals to be measured, stored in the memory 13, by which is 
obtained a complex base band signal l S( Qs having errors eliminated therefrom and processed by the receiving filter. 
That is, the signal by which the signal read out of the memory 13 is multiplied is a signal obtained by subtracting the 
carrier frequency error from the frequency of a locally generated intermediate-frequency signal and correcting its initial 
phasa In other words, this signal is exp(-j(co-co 1 -a> 2 )-9 (where 9 is the phase corresponding tox^T 2 +8 o )• ln this 

40 instance, the initial phase is the sum of the phase corresponding to the time that is the sum of a deviation x<\ from the 
symbol decision point and a deviation x 2 from the pilot PN pattern and an error of the initial phase 8 0 by the carrier fre- 
quency error. The low-pass filter for transforming the signal read out of the memory 13 to the base band signal has the 
same characteristic as that of the receiving filter (a complementary filter). The data to be measured starting at the 
beginning of the PN code has a length corresponding to 64 by 20 chips in this example. 

45 Then, an error vector is calculated from the complex base band signal l 5 ,Q 5 thus orthogonal-transformed so as to 
correct the frequency error and so forth and the error vector is used to calculate modulation accuracy, amplitude accu- 
racy, phase error and IQ origin offset. Sufficient accuracy for the calculation of the waveform quality p can be obtained 
with the base band signal l 5 ,Q 5 but amplitude and initial phase errors are important for the calculation of the modulation 
accuracy. Hence, the amplitude error, the initial phase and the carrier frequency error are estimated from the complex 

so base band signal l 5 .Q 5 in a manner to meet a modulation accuracy computing standard, and these estimated parame- 
ters are used to correct the complex base band signal l 5 ,Q5, followed by the calculation of the modulation accuracy and 
other parameters. 

That is to say, these parameters are estimated for the base band signal l 5 ,Q s in a frequency error, phase and ampli- 
tude estimating section 24. Letting the number of the signal data sequence be represented by k. the amplitude of the 
55 ideal signal R|A by I. the amplitude and phase components of the base band signal l 5 ,Q 5 by a(k) and 6(k) and the 
phase component and the amplitude of the ideal signal R If Q| by 9 R (k) and a R (k), a carrier frequency error Ac* the initial 
phase 6 0 and an amplitude coefficient b 0 are calculated by the following equations. 
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Aco = {(Za(k) • [6(k) - e R (k)]) . (Za(k) • k) - (Za(k)) • (Za(k) • [e(k) - 
e R (k)] - k)}/{(2a(k)) 2 - (Za(k)) • (SaW • k 2 )} 

s e 0 = {(za(k) • k) • (za(k) • [e(k) - e R (k)] • k) - (Za(k) • (e(k) - 

e R (k}]) • (za(k)}/{(za(k) • k) 2 - (za(k)) • (za(k) • k 2 )} 

b o = {Za(k) • a R (k) * cos[6 R (k) - 6(k) + Aa> • k + 6 0 ]}/Za(k) 2 

where Z indicates the summation from 1 of k to the last k (=N) of the input signal. 

These calculated values At* 0 0 and b 0 are used to correct the input complex base band signal l 5 ,Qs by the following 
equation in a frequency, phase amplitude correcting section 25. 

75 

1 6 + jQ 6 = b 0 (1 5 + jQ s )exp(-j( A© • k + e 0 ) 

The thus corrected base band signal fe.Qe and the ideal signal R|.Q| are used to conduct the following equation in 
an error vector and waveform quality p calculating section 26. 
20 In the first place, squares of the amplitude of the error vector at symbol points (every chip period) are all added 
together by the following equation. 

Z|e(k)| 2 = Z[(l 6 (k) -R,(k)) 2 +(Q e (k). R Q (k)) 2 ] 

25 where Z indcates the summation from k=1 to N. 

This value is divided by the number N of additions and the square root of the divided value is computed, followed 
by calculating the following equation to express it as a percentage with respect to the amplitude of the ideal signal. 

s r a100x V{Z|e(k)| 2 /N)%rms 

30 

Next to obtain the IQ origin offset added values of the error vector in the I- and Q-axis directions are calculated by 
the following equations. 

Xe.M-ZCIM-R.M) 

35 

£e Q (k) = X(l(k)-R Q (k)) 

where Z indcates the summation from k=l to N. 

The added value in each axis is divided by the number N of additions to obtain the mean value, which is converted 
40 to a dB value. Since the amplitude of the ideal signal is 1, letting the amplitude be represented by 0 dB, the IQ origin 
offset is given by the following equation. 

O 0 = 20 x log 10 {[V(Z*,(k)) 2 + (Z* Q (k)) 2 ]/N 2 } 

45 Next, the amplitude error is calculated. This begins with subtracting the amplitude 1 of the ideal signal from the 
amplitude of the measured signal l 6 ,Q 6 , followed by squaring the remainder. 

Z e m (k) 2 = Z(V(I S 2 (k) +Q 6 2 (k)) - 1) 2 

so where Z indicates the summation from k= 1 to N. This is converted by the following equation to the %rms unit to compute 
the amplitude error. 

emr = 100x V(Ze m 2 (k)/N) 

55 After this, the phase error is calculated. In this instance, the phase of the ideal signal Rj.Rq is subtracted from the 
phase of the measured signal l 6 ,Q 6 and the remainder is squared. In this subtraction of phase, the phase difference 
may sometimes be calculated in the opposite direction when the phase is near the 180° axis; so that the phase differ- 
ence is corrected when it is larger than 180°. 
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Xe p 2 (k) = z(arctan(Q(k)/l(k) - e R (k)) 2 



where X indicates the summation from k=1 to N. This is converted by the following equation to a unit deg rms. 



5 



Cpr »(180/*)x V(Ze p 2 (k)/N) 



The waveform quality p is calculated by the following equations based on its definition. 



ReQ) = S(l e (k) • R R W + Q 8 (k) • R Q (k)) 



10 



lm(j) => Z(Q 6 (k) - R R (k) - 1 6 (k) • R Q (k}) 



where Z indicates the summation from k»1+64(j-1) to 64+64Q-1). 



15 



P = S(ReO) • ReG) + lm(j) • lm(j)/[64 • S(l 3 2 (k) + Q 8 2 (k))] 



where 2 in the denominator indicates the summation from k=1 to 64N and X in the numerator the summation from k=1 
toN. 

The frequency error <o^ estimated in the frequency estimating section 15 ( the frequency error in the frequency 
20 error and phase estimating section 19 and the frequency error A© in the frequency error, phase and amplitude estimat- 
ing section 24 are added together in a common unit in a frequency error calculating section 27 to obtain a frequency 
error 004. 

From the PN synchronized ideal signal generating section 1 8 a difference t 1 between the trigger from the base sta- 
tion and the beginning of the PN code sequence of the received input signal is obtained, which is provided to a time 

25 alignment error calculating section 28 for addition to a signal delay x 2 estimated in the signal delay estimating section 
22. The difference (1/2 chip) between reference timing defined by the standard and the symbol decision point and the 
pre-measured amount of delay by hardware in the measuring system are used to correct the added value T1+T2 t0 
obtain a time alignment error x 3 . 

The signal l 6 ,Q 6 and the parameters p, ^ ep P O 0i oo 4 and x 3 obtained as described above are displayed on a 

30 • display 29. In an error and ideal signal estimating section 30 made up of the orthogonal transform section 14 through 
the signal delay estimating section 22 in Fig. 2, a frequency error and a timing error of the input signal and an ideal ref- 
erence signal are estimated from the input signal whose length is significantly shorter than that of the signal to be meas- 
ured, for example, about 400 chips (symbols), and the estimated values are used to process the entire input signal to 
be measured in the processing sections 23 through 26. 

35 While in the above intersymbol interference is eliminated from the input signal in the PN synchronized ideal signal 
generating section 18, a signal with intersymbol interference or data other than the symbol decision point may some- 
times be needed in the case where the modulation system used is not that of the CDMA digital cellular telephone. In 
such an instance, PN patterns processed by the sending filter are held in a memory Although the present invention has 
been described as being applied to the measurement of the waveform quality of the modulated signal of the spread 

40 spectrum system, the invention is also applicable to the measurement of the waveform quality of modulated signals of 
other systems-this can be achieved usually by utilizing a predetermined known word portion such as a synchronization 
word. In the case of using part of the PN code as transmitted data, the known data can be used in synchronization with 
the PN code portion. The receiving filter 1 7 in Fig. 2 may sometimes be omitted. 

As described above, according to the present invention, since the ideal reference signal is generated using a known 

45 portion of the symbol in the input signal, not from input signal demodulated data, no measurement error will not be 
caused by a demodulation error. In this respect the present invention permits effective measurement of modulation 
parameters of signals prone to introduce demodulation errors. 

Further, according to the present invention, the frequency error is relatively easily obtained using a portion of the 
input signal, for example, about 400 symbols and the frequency error is used to correct the input signal. The ideal ref- 

50 erence signal is easily derived from the corrected input signal I2.Q2 and the ideal reference signal thus obtained and 
the corrected signal l 2 ,Q2 are used to obtain frequency and timing (phase) errors relatively easily. Thereafter, the input 
signal is subjected to orthogonal transformation to obtain a complex base band signal le,Qe so that these errors are 
removed, and modulation parameters are derived from the complex base band signal and the ideal reference signal 
R|,Rq. The computational quantity in the present invention is smaller than in the case of increasing the measurement 

55 accuracy by repeating the production of all modulation parameters as in the afore 

mentioned U.S patent. Repeating the calculation of processed data, there is a fear that calculation errors accumu- 
late, making the calculated value rather inaccurate. According to the present invention, however, the input signal is sub- 
jected to orthogonal transformation in the processing section 23 and its output is processed only once; hence, there is 
no accumulation of calculation errors by repeating calculations. 
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By obtaining the frequency error, timing error and ideal reference signal with a signal shorter than the input signal 
to be measured, the computation complexity can be further reduced. 

In the case of performing the filter processing, the operation therefor is particularly time-consuming, but the com- 
putational complexity needed in the present invention can be reduced far smaller than in the prior art which always uses 
5 the entire signal to be measured. 

In the embodiment described above, since the frequency error is corrected prior to the receiving filter processing, 
part of the signal is not cut off by the receiving filter processing, ensuring accurate measurement 

Since the CDMA mobile communication signal is spread by the PN code, if each chip is used as the symbol to be 
measured in the conventional method, the number of symbols becomes, for example, six-fold except the pilot signal and 
10 the computing time for measurement also becomes six-fold. However, the present invention is far smaller in computa- 
tion complexity as described previously, and hence it is advantageous over the prior art. 

According to the present invention, the input n-phase PSK signal is converted to a complex base band signal and 
it is raised to the n-th power as shown in Fig. 3. Accordingly, the carrier frequency of the complex base band signal is 
an error frequency of the input n-phase PSK signal, and this is a significantly small value and the sampling frequency 
is of the signal raised to the n-th power satisfies the concfition of the sampling theorem without fail, making it possible to 
obtain the carrier frequency with accuracy. 

In the synchronization of the spreading code with the input signal, the symbol decision point is estimated and the 
correlation is calculated for only samples near the symbol decision point as described previously with reference to Figs. 
4 and 5, the computational complexity required is smaller than in the case of calculating the correlation for every sam- 
20 pie. In the case of the pilot signal, in particular, since the correlation with the samples near the symbol decision points 
is calculated for only the real or imaginary part of the input complex base band signal, the computational complexity is 
sharply reduced as compared with that in the prior art 

in the operation in the signal delay estimating section 22, that is, in the estimation of the symbol timing, the products 
A, B and C of the value of cross correlation between the input complex base band signal Z(k) and the ideal signal R(k) 
25 and the filter coefficients a m , b m and c m are computed as indicated by Eqs. (16), (17) and (18) and x is calculated by 
Eq. (20). Accordingly, the value x can be obtained with high accuracy without involving any interpolating operation and 
hence in a short time (with small computational complexity). In contrast to this, according to the conventional method of 
Fig. 1 . the calculation of the value x with high resolution and consequently with high accuracy, requires the reduction of 
the sampling interval by performing an interpolation operation. That is, the conventional method repeats the calculation 
30 while shifting the input signal and the ideal signal relative to each other, and hence it consumes much time for the cal- 
culation. 

It will be apparent that many modifications and variations may be effected without departing from the scope of the 
novel concepts of the present invention. 

35 Claims 

1. A modulation parameter measuring method, in which: a digital sequence of an input digital quadrature-modulated 
signal stored in a memory is read out therefrom and subjected to orthogonal transformation into a complex base 
band signal; an ideal reference signal of said complex base band signal, synchronized with its symbols, is gener- 
40 ated; a frequency error, an initial phase error and a timing error of said complex base band signal are estimated; 
said frequency error, said initial phase error and said timing error are used to correct said complex base band sig- 
nal; and said corrected complex base band signal and said ideal reference signal are used to measure modulation 
parameters of said input digital quadrature-modulated signal; 
CHARACTERIZED BY: 

A first step of reading out part of said digital sequence from said memory and subjecting it to orthogonal trans- 
formation into a first complex base band signal; 

a second step (30) of estimating a frequency error (hereinafter referred to as a first frequency error), an initial 
phase error (hereinafter referred to as a first initial phase error) and a timing error (hereinafter referred to as a 
so first timing error) of said first complex base band signal and an ideal reference signal; and 

a third step (2!3) of obtaining said corrected complex base band signal by orthogonal transformation of the 
entire digital sequence necessary for measurement stored in said memory in such a manner as to correct said 
first frequency error, said first initial phase error and said first timing error. 

55 2. The method of claim 1 . wherein said second step comprises: a seventh step (1 8) of generating said ideal reference 
signal by using said first complex base band signal; and 

an eighth step (19) of estimating a frequency error hereinafter referred to as a second frequency error) and 
initial phase error (hereinafter referred to as a second initial phase error) from said first complex base band signal 
and said ideal reference signal, said second frequency error and said second initial phase error being respectively 
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said first frequency error and said first initial phase error in said third step. 

3. The method of claim 2, wherein said eighth step is a step of calculating said first frequency error by calculating the 
rate of change of the phase difference between said ideal reference signal and said first complex base band signal 

5 with respect to the lapse of time, an initial value of the phase difference between said ideal reference signal and 

said first complex base band signal being said first initial phase error. 

4. The method of claim 2, wherein said seventh step comprises: a symbol decision point estimating step (131) of esti- 
mating symbol decision points from said first complex base band signal oversampled with respect to a symbol 

10 period; a decimating step (1 37) of decimating samples close to said estimated symbol decision points from the one 
of two components of said first complex base band signal; a correlation calculating step (138, 139) of calculating 
the correlation between said decimated samples and a complex local signal; and a step of obtaining said ideal ref- 
erence signal by controlling the delay of said complex local signal so that the value of said correlation is maximized. 

is 5. The method of claim 4, wherein said symbol decision point estimating step (131): a step (132, 133, 134) of calcu- 
lating the magnitude of each sample of said first complex base band signal; a step (1 35) of performing discrete Fou- 
rier transform processing of a sequence of said magnitudes; and a step (136) of calculating said symbol decision 
points from the phase of a symbol period component contained in the result of said discrete Fourier transform 
processing. 

20 

6, The method of claim 5, wherein said discrete Fourier transform processing is performed by multiplying said magni- 
tude sequence by instantaneous values of sine and cosine waveforms for each sample period. 

7. The method of claim 2, wherein said second step comprises a ninth step (1 5) of estimating a frequency error (here- 
25 inafter referred to as a third frequency error) of said first complex base band signal, and a tenth step (1 6) of correct- 
ing said first complex base band signal with said third frequency error to obtain said first complex base band signal 
for use in said seventh step, the sum of said third frequency error and said second frequency error obtained in said 
eighth step being said first frequency error in said third step. 

30 8. The method of claim 7, wherein said ninth step (15) comprises: a step (124) of raising said first complex base band 
signal of n-phase modulation to the n-th power; a step (125) of performing fast Fourier transform processing of said 
first complex base band signal raised to the n-th power; and a step (126) of frequency-cfividing a frequency com- 
ponent in the output of said fast Fourier transform processing, which provides the maximum peak, down to 1/m to 
obtain said third frequency error. 

35 

9. The method of claim 2, which further comprises: an eleventh step (21) of correcting said first complex base band 
signal with said second frequency error and said second initial phase error, and a twelfth step (22) of calculating a 
timing error (hereinafter referred to as a third timing error) between said first complex base band signal corrected 
in said eleventh step and said ideal reference signal, said second timing error being said first timing error in said 

40 third step. 

10. The method of claim 9, wherein said twelfth step comprises: a step (151) of calculating the cross correlation 
between said first complex base band signal corrected in said eleventh step and said ideal reference signal; a step 
(152, 153, 154) of performing filter processing of said cross correlation by three Nyquist filter characteristics each 

45 approximated as a function of said third timing error; and a step (155) of calculating said third timing error from the 
results of filter processing using said three Nyquist filter characteristics. 

11. The method of any one of claims 1, 2, 3, 4 f 5, 6, 9 and 10. wherein said second step further comprises: a thirteen 
step (1 5) of a frequency error (hereinafter referred to as a third frequency error) of said first complex base band sig- 

so nal; a fourteenth step (16) of correcting said first complex base band signal with said third frequency error; and a 
fifteenth step (1 7) of subjecting said first complex base band signal corrected in said fourteenth step to processing 
a filter of a predetermined characteristic to obtain said first complex base band signal for use in said second step. 

12. The method of any one of claims 1 through 11, which further comprises a step of detecting the time difference 
55 between a trigger from the outside and said ideal reference signal as a time alignment error, the sum of said first 

timing error in said second step and said time alignment error being said first timing error in said third step. 

13. The method of claim 12, which further comprises the steps of: reading said input digital quadrature-modulated sig- 
nal into a memory using said trigger as the reference; and detecting, as said time alignment error, the difference 
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between the address of said memory corresponding to said trigger and the address corresponding to the beginning 
of said ideal reference signal. 

14. The method of claim 1 4, wherein part of said digital sequence is about 400 symbols long. 

15. The method of any one of claims 1 through 14, wherein modulation data of said input digital quadrature-modulated 
signal is preknown. 

16. The method of any one of claims 1 through 15, wherein said input digital quadrature-modulated signal is a trans- 
mitted signal modulated by a pilot PN code in a CDMA mobile communication system. 
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(57) The signal to be measured is converted into 
digital form and read into a memory 13 using a trigger 
as the reference, and a signal of a length 1/5 that of the 
signal to be measured, in the vicinity of the trigger, is 
orthogonal-transformed (14). A frequency error co 1 of 
the transformed output l 1t Q-| is estimated (15) and is 
used to correct the output l|,Qi (16). The corrected out- 
put I2.Q2 is subjected to processing by a receiving filter 
(17), and the filtered output 1 3 ,Q 3 is synchronized with a 
PfsJ code to obtain an ideal signal R|.Rq and a phase dif- 
ference -c 1 is calculated. The signal l 3 ,Q 3 and the ideal 
signal R|,Rq are used to obtain a frequency error g>2 
and a phase error T2 of the signal l3,Q 3 . All signals to be 
measured, stored in the memory 13, are subjected to 
orthogonal transformation (23) so that the frequency 
error m^o^ and the phase error x 1 +x 2 are removed, 
and the transformed output l 5 ,Q 5 and the signal R|,Rq 
are used to calculate modulation parameters as is the 
case with the prior art. 
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